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ABSTRACT 

The primordial power spectrum describes the initial perturbations in the Universe which eventually 
grew into the large-scale structure we observe today, and thereby provides an indirect probe of infla- 
tion or other structure-formation mechanisms. In this paper we will investigate the best scales the 
primordial power spectrum can be probed, in accordance with the knowledge about other cosmological 
parameters such as f^, f2 c , J7a, h and t. The aim is to find the most informative way of measuring 
the primordial power spectrum at different length scales, using different types of surveys and the 
information they provide for the desired cosmological parameters. We will find the optimal binning of 
the primordial power spectrum for this purpose, by making use of the Fisher matrix formalism. We 
will then find a statistically orthogonal basis for a set of cosmological parameters, mentioned above, 
and a set of bins of the primordial power spectrum to investigate the correlation between the two 
sets. For this purpose we make use of principal component analysis and Hermitian square root of the 
Fisher matrix. The surveys used in this project are Planck and SDSS(BRG), but the formalism can 
easily be extended to any windowed measurements of the perturbation spectrum. 
Subject headings: Cosmic microwave background — cosmological parameters — early universe - 
large-scale structure of universe 



1. INTRODUCTION 

The primordial Power Spectrum (PS) probes the 
physics of structure formation in the early Universe. In 
particular, inflation provides a paradigm for early Uni- 
verse physics in accord with current cosmological obser- 
vations. Simple models of inflation predict an almost 
Gaussian distribution of adiabatic perturbations with a 
scale-invariant spectrum (i.e., P(k) cx k). However, there 
are other possibilities: for example, there could be more 
than one scalar field during inflation and this would pre- 
dict a different spectrum and possibly a different dis- 
tribution of fluctuations. The details of inflation are 
presently unknown to us. Therefore, determining the 
primordial PS would give us better intuition about the 
early Universe. 

We use different surveys — of the Cosmic Microwave 
Background, of the galaxy power spectrum, of velocity 
fields, etc. — to constrain the primordial PS. However, 
the spectra of these surveys are jointly sensitive to cosmo- 
logical parameters (which, we will collectively call 0j)and 
the primordial PS. Hence there is a statistical degeneracy 
between the two. The aim here is to explore, as the data 
improve, what new information can be learnt about the 
primordial PS and what exactly needs to be improved to 
better constrain the primordial PS. The motivation is to 
test the assumptions about the initial conditions besides 
getting better constraints on parameters based on the 
same set of assumptions. Therefore, knowing the degen- 
eracy between the cosmological parameters and primor- 
dial PS, we want to investigate the scales the primordial 
PS can be probed best with future experiments. 

The outcome of different surveys is usually a type of 
PS that is a convolution of the primordial PS, whatever 
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form it may have, and (the square of) the transfer func- 
tion of the particular type of survey, which holds the 
cosmological parameters. Here I list some examples; 

• For galaxy surveys, the PS is related to the primor- 
dial PS through the matter PS, P$(k), as 

P g (k) = b 2 (k)P s (k) ~ b 2 (k)2ir 2 kT 2 (k)A 2 c (k) , (I) 

where is the primordial PS and T(k) is the mat- 
ter transfer function and b(k) is the bias. 

• For CMB surveys, the angular PS is 

C e = 4ir d\nkA 2 e (k)A 2 c (k) , (2) 
Jo 

where £ is the related to the angular scale on sky 
via £ ~ 180° /0 and Af(fc) is the angular transfer 
function of the radiation anisotropics. 

Here, we define the primoridial curvature power spec- 
trum, which we parameterize as A 2 (k) = i4.(fc/0.05) n * -1 . 
A is the amplitude and n s is the spectral index. The 
notation refers to the g auge-invariant curvature pertur- 
bation C (|Bardeerill98oh . 

Other types of power spectra, such as the weak lensing 
and peculiar velocity power spectra, have similar forms; 
they depend on the cosmological parameters, through a 
transfer function, and the primordial PS. These different 
power spectra probe different scales with different accu- 
racies. One survey can, therefore, help fill the gaps in 
other surveys and all together they are expected to im- 
prove information, especially on the overlapping scales. 
This means combining surveys can help us choose nar- 
rower bins and hence investigate the primordial PS to a 
greater resolution. 
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One common method for error estimation is to use a 
Fisher matrix analysis. The Fisher matrix is generally 
used to determine the sensitivity of a particular survey 
to a set of parameters and has been largely used for fore- 
casting and optimisation. The Fisher matrix is the en- 
semble average of the curvature of a function T (i.e. it 
is the average of the curvature over many realisations of 
signal and noise) 
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The Fisher matrix allows us to estimate the errors on 
parameters without having to cover the whole parameter 
space. Hence, the inverse of the Fisher matrix is a crude 
estimate of covariance matrix of the parameters, by anal- 
ogy with a Gaussian distribution in the 9 a , for whic h this 
would be exact. The authors of dBond et al.lll998f ) have 
compared the Fisher matrix analysis with the full likeli- 
hood function analysis and found there was great agree- 
ment between the two methods if the likelihood function 
is approximately Gaussian near the peak. The Cramer- 
Rao inequality states that the smallest non-marginalised 
error measured for the parameters by any unbiased esti- 
mator (such as the maximum likelihood) is 1 / \[F 1 . The 

marginalised 2 one-sigma error is \J (F^ 1 )^ for parame- 
ter a. 

The Fisher matrix for CMB surveys is given by 
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where Ct is the angular PS, w is the weight defined as 
(AQ<7„) with Af2 being the real space pixel size and 

ct 2 , being the noise per pixel, e~ e a is the window func- 



tion 3 for a Gaussian beam (where a = 9f w hm /V8 In 2) 
and f s ky is the fraction of the sky observed. The fac- 
tor fskyi^i + 1) gives the number of independent modes 
at a given wavenumber; the term proportional to Ci is 
the sample (or cosmic) variance contribution, and the 

w~ 1 e t a term is the noise contribution. Note that the 
diagonal form for the matrix implies diagonal (uncorre- 
cted) errors on the CiS. To find errors on other param- 
eters, we use the Jacobian 
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where 9 a and Op are different parameters. 

For a volume -limited galaxy survey the Fisher matrix 
(|Tegmarklll997ft is 
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(27r) 2 (P„ + l/n) 2 



(6) 



where V is the total volume of the survey, n is the num- 
ber density of the survey (N to t/V), P n is the galaxy PS 

It should be noted that the Cramer-Rao inequality is a state- 
ment about the so-called "Frequentist" confidence intervals and is 
not strictly applicable to "Bayesian" errors. 

2 Integration of the joint probability over other parameters. 

3 This damps power on larger is; as we get closer to the resolution 
limit of the survey CfS start to correlate. 

4 Note that this equation only applies to linear regime, as non- 
linearities impose non-Gaussianities. 



in each k n bin and Afc is the binwidth. Similar to the 
CMB power spectrum case, fc 2 AkV counts the number of 
modes, P n gives the sample variance, and 1/n the noise 
variance due to Poisson counting errors. This, again, 
gives us the errors on the galaxy PS and we use the Ja- 
cobian to get the errors on other parameters 



F a f3 — J]] F n 



dP(k n ) dP(k n 
d6 a Wr 



(7) 



Fisher matrices for different surveys can easily be com- 
bined by a simple summation F = Fg a i ax y + Fpjyjg. 

This is because they are proportional to the log of the 
likelihood function and we multiply likelihoods to com- 
bine them. Equivalently, we can think of them as the 
weights (inverse noise variance) of the experiments, which 
add for a Gaussian distribution. The nonzero correlation 
between the parameters in the covariance matrix makes 
interpreting the errors somewhat more difficult than the 
uncorrelated case. We will discuss various methods for 
decorrelating the power spectra and cosmological param- 
eters. 

2. METHOD 

The aim is to investigate the primordial PS in a "non- 
parametric" way (we use quotations remarks to remind 
the reader that "non-parametric" merely means that we 
use a very general model, potentially with a very large 
number of parameters). Therefore, we assume a top- hat 
binning of the primordial PS 



A 2 (fc) =J2 w B(k)Q B 



(8) 



where Qb is the power in each bin B and wb = 1 if 
fc G B and otherwise. The cosmological parameters 
under investigation are (and of the form) O c , f2f,, f^A, 
h, t and n s . The reason for inclusion of n s in the pa- 
rameters is to allow for a consistency check. In this set- 
ting we do not expect to see any correlation between n s 
and the bins of primordial PS. Inclusion of n s in the pa- 
rameter space only makes minute changes to our results 
and can be ignored. We will choose a ge ometrically flat 
(adia batic) ACDM model with WMAP5 (jDunklev et all 
12008ft values for the parameters; n s = 0.963±0.0145(with 
zero running), fi m = 0.214 ± 0.027, Q b = 0.044 ± 0.003, 
0\ = 0.742 ± 0.03, t = 0.087 ± 0.017 and h = 0.719 ± 
0.0265, where H = lOOMtm^Mpc -1 . Q. v = 0.0 was 
chosen, as massive neutrinos in troduce some difficul- 
ties in the Fisher matrix analysis (jEisc nstein ct al.lll999T ) 
and therefore were ignored for now. CMBfast software 
(|Seliak fc Zaldarriaeal H99l was used for the calcula- 
tions. The surveys chosen for this initial investigation are 
the projected results from the SDSS Bright Red Galax- 
ics(BRG) 5 sample and the Planck Surveyor CMB Power 
Spectrum 6 . 

3. GALAXY SURVEYS — SDSS(BRG) 

5 These are bright galaxies, which means the survey can be quite 
deep, with z ~ 0.25 — 0.5. Also, these trace the elliptical galaxies, 
which are thought to be better tracers of mass at this rcdshift 
range 

6 
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A galaxy PS is related to the matter PS via a 
parameter called bias — equation [1] For the BRG 
sample of SDSS, this is assumed linear and scale- 
independent with the form P g = b 2 P$, where & is the 
bias and approximately equal to 2.0 (IMann et alj|1998t 
Scherrer fc W^einberdll998t lHutsil[200alSeliak fc Warrenl 



the primordial PS bins, the derivative becomes 7 



2004) . The survey specifications fo r BRG sample are 
n = 10 5 /V and V = (l/i^Gpc) 3 (IGunn fc Weinberg! 
[19951) . 

For the 9i the derivatives in the Jacobian were obtained 
numerically using the Taylor expansion 



PWi) = p(e ) + (^)A(<?o 



(9) 



The width and direction of the step are quite important 
here. A two-sided derivative was chosen, so that the 
derivative is centred on the default value 6q, with a step 
size of A(0j)/2 on each side. This is accurate to 2 nd order 
in A(0j) (a one-sided derivative would be at a slightly 
shifted place of Oj + A(6j)/2, and is only accurate to I s * 
order (|Eisenstein et al.l Il999f )). The width of the step 
should be small enough to give accurate results and yet 
big enough to avoid numerical difficulties. This was taken 
to be a 5% variation, therefore a 2.5% width on each 
side. Other studies have shown that this turns out to 
be the best step size, giving the most accurate results 
(|Eisenstein et al.lll999ft . 

For the primordial PS bins, the derivative is propor- 
tional to the matter transfer function 



dPgjk) 
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(10) 



where fc and k! refer to the bins. The fc-range for SDSS 
is 0.006 < fc/^Mpc -1 ) < 0.1. The minimum value 
for the wavenumber, k m i n , is obtained from the largest 
scale of the survey — (2-rr/V 1 ^ 3 ). Its maximum value, 
kmax, is chosen to avoid non-linearities. Sim ulations of 
a very similar flat model (|Meiksin et al.l[l998l ) suggested 
a k max of O.l/iMpc -1 . This is also very close to the 
scale at which departure s from linear theory was seen by 
IPercival fc White! ((20081) . 

The derivatives in the Jacobian need to be averaged 
into bins. Later we will explain the criteria for choosing 
the widths and locations of the bins. 

4. CMB SURVEYS — PLANCK 

One thing to note in this case is that the output of 
CMBfast is of the form C e = [£{£ + 1)/2tt]CV, so the 
CMB Fisher matrix, equation is multiplied by this 
factor. The specifications for Planck HFI (y = 100GHz) 
are 9 fwhm = 10.7' = 0.0 03115 radians, a pix = 1.7 x 10~ 6 , 
to" 1 = 0.028 x 10" 15 (jDelabrouille et all [1991 . The 
derivatives in the Jacobian were again obtained numeri- 
cally by Taylor expansion 



c l {e i )=c l {e ) + {^-)A{e i ) 



(11) 



The same arguments as in the SDSS case applies for 
the width and direction of the step here. In the case of 
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This needs to be averaged into k bins, which will be 
explained later on. The chosen fc-range for Planck is 
0.0001 < k/(hMpc^) < 0.1, where k min was obtained 
from k m in = (-min/dA = 2/4a, where c?^ is the angular 
diameter distance to t he surface of last sc attering ob- 
tained to be ~ 14Gpc (jDunklev et al.ll2008D . 

5. SDSS & PLANCK 

As explained above, to combine data from different sur- 
veys, we can add the Fisher matrices obtained for each 
of them. We expect to see an improvement on the errors 
of both the bins and cosmological parameters. Equiva- 
lently, this will enable us to have narrower bins without 
sacrificing Signal-to-Noise per bin. 

6. OPTIMAL BINNING 

As explained before, a set of primordial PS bins are 
part of our parameter space. In this section we will ex- 
plain how these bins are chosen. For our purposes the 
bins need to have the same amount of contribution to the 
Fisher matrix which means they need to have the same 
S/N. We take the signal in each bin to be the amplitude 
of the primordial PS in that bin and the noise to be given 
by the inverse of the square root of the diagonal elements 
of the Fisher matrix. For this, we construct a signal vec- 
tor, S, which contains the amplitude of the primordial 
PS for all the bins and the values of the cosmological 
parameters. Wc weight our Fisher matrix by this vector 



(13) 

where there is no Einstein summation. This now gives 
us a (S/N) 2 matrix, where the square root of its diago- 
nal elements are the S/N for the bins, and the weighted 
errors for #;S. It is worth to emphasise that it is this 
(S/N) 2 Fisher matrix that will be diagonalised later on. 

For the SDSS case, we start by having the maximum 
number of bins possible in our fc-range. The usual prop- 
erties of the Fourier transform imply that the scale of the 
survey not only determines k m i n , but also puts a limit on 
our resolution: k m i n — (Ak) m i n = (2n/V 1 ^ 3 ); narrower 
bins would become highly correlated. Therefore, we set 
up a series of bins with this minimum binwidth in our 
fc-range. We then construct a Fisher matrix for this set 
of bins (and 9is) and weight it by the signal vector, S, 
for this set. With this binning adopted, the S/N values 
range from 3.7 in the first bin to 35.1 in the last bin. 
Knowing that the binwidths chosen are the minimum 
possible and that increasing binwidths will increase the 
S/N value, we conclude that the bin with the maximum 
S/N cannot be changed and hence we make other bins 
wide enough to reach the S/N in this bin. To obtain this 
optimal binning we start an iteration; smaller bins of size 
~ (Afc) m i„/6 are combined until the S/N are equal to 
15% of the maximum S/N: 



Max(S/A0 - (S/N), 



0.15. 



(14) 



Max(S/iV) 

7 To obtain A((k), CMBfast needed to be altered to give the 
radiation transfer functions at all is. Then, for each £, this was 
interpolated in k. 



4 



where i refers to the bins. This gives us 8 bins with their 
S/N ranging 30 - 35. 

For Planck, the bins are obtained so that their S/N 
matches that of SDSS. The reason for applying this cri- 
teria to Planck is to allow for a fair comparison between 
the results from SDSS and Planck. This criteria gives us 
a total of 23 bins for Planck. 

In the case of the combined Planck and SDSS we re- 
quire only that the S/N of the bins are equal to 50%. This 
now gives us the optimal resolution of the primordial PS 
we can achieve from SDSS and Planck. We have a total 
of 48 bins with S/N being in the vicinity of ~ 20.0 and, 
therefore, still comparable to the S/N values in the other 
cases. 

It is worth reminding the reader that an aternative 
way to determine the binning would be to take the 
marginalised errors as the noise. This would be obtained 
by inverting the Fisher matrix in each iteration loop to 
get the covariance matrix, which gives the marginalised 
variances of the bins and 0jS. We would then take the 
sub-block of this covariance matrix that refers to the bins 
only, and invert it to get a marginalised Fisher matrix 
for the bins. We would then feed this Fisher matrix into 
equation 1131 However, this method could not be imple- 
mented because the SDSS Fisher matrix is not invert- 
ible; the SDSS Fisher matrix is not a positive definite 
matrix because it is asked to estimate too many param- 
eters. There are a total of n data points (n galaxy PS 
bins) and we are asking these to predict n + m parame- 
ters (n primordial PS bins and m #js). Also, note that 
whichever of the methods presented uses the correlated 
errors as the noise. We now discuss the decorrelation of 
the parameters. 

7. DECORRELATING THE PARAMETERS 

7.1. Principal Component Analysis 

One popular method to overcome the correlation be- 
tween the parameters is to perform principal component 
analysis (PCA); the covariance matrix is a symmetric 
n x n matrix and therefore, can be diagonalised using 
its eigenvectors. This has the form C = E T AE, where 
C is the covariance matrix, E is an orthogonal matrix 
with the eigenvectors of C as its rows and A is the di- 
agonal matrix with the eigenvalues of C as its diagonal 
elements 8 . This constructs a new set of variables X that 
are orthogonal to each other and are a linear combination 
of the old parameters O, through the eigenvectors 



X = EO 



(15) 



The Xi are called the principal components of the ex- 
periment and are ordered so that X\ and X n are the 
best and worst measured components respectively In 
this construction, the eigenvalues are the variances of the 
new parameters so that X\ has the smallest eigenvalue 
and X n has the largest one. The eigenvectors have been 
normalised so that 53 . e| = 1, where ejS are the elements 
of Bj. We list some properties of PCA below; 

• The main point of PCA is to assess the degenera- 
cies amongst the parameters that are not resolved 

8 It is common to construct the covariance matrix for the PCA. 
However, Fisher matrix can be used instead; eigenvectors stay the 
same, but eigenvalues are inversed. 



by the experiments, be they fundamental like the 
cosmic variance or due to the noise and coverage 
of the experiment. In our case, it will especially 
help us to see the correlation amongst the bins of 
the primordial PS, and between the bins and the 
cosmological parameters. 

• The eigenvalues obtained measure the performance 
of the experiment — a larger number of small eigen- 
values means a better experiment. Another mea- 
sure of the performance of the experiments is to 
see how they mix physically independent parame- 
ters such as, say, n s , the spectral index, and f2&. 
This sort of mixture may be improved by improv- 
ing the experiment. 9 

All the above points may be summed up to conclude 
that in a perfect setting we would expect a one-to-one 
relation between the old and the new parameters. This 
means that we would see only one of the old parameters 
to contribute completely to one of the new parameters, 
with zero contribution from the other old parameters. 

Note that the principal components obtained are not 
unique and depend on the form of the variables (e.g., 
whether we use f2t or log fib), as well as where they are 
evaluated. 

7.2. Hermitian Square Root 

Another approach to remove the correlations between 
the uncertainties is to use the Hermitian square root 
of the Fisher ma trix as a linear transformation on the 
parameter space (jBond et all 119981 ; [Hamilton! ll997alTbT ). 
This transformation matrix is obtained by 



F 1 / 2 = E T A 1 / 2 E . 



(16) 



where, like before, E is the eigenvector matrix and A 
is a diagonal matrix containing the eigenvalues. It has 
the property F = F 1/2 F 1/2 = (F 1/2 ) T F 1/2 and therefore 
the condition (F- 1 / 2 )F(F- 1 / 2 ) = (F-i/3)T F ( F -i/a) = 
diag is satisfied. Unlike PCA, F 1 / 2 docs not give us an 
orthogonal basis and instead, it can be thought of as 
'window functions' for the primordial PS. We define a 
window matrix by 



(JTl/2) 



(17) 



which satisfies the normalisation condition J~] H nrn = 1 . 
Hence the windowed PS is defined as 



(18) 



where P(k n ) is the original primordial PS. Note that this 
windowed PS is not a physically motivated PS and it 
is just constructed for a visual presentation and under- 
standing of the underlying correlations. However, in a 
perfect setting we would expect this windowed PS to be 

9 However, the so-called 'geometr i cal degeneracy' 
flZaldarriaea et al.l 119971; lEfstathiou fc Bond! IT999D will not 
be improved by improving the experiments; two models with same 
primordial PS, the same matter content, and the same comoving 
distance to the surface of last scattering produce identical CMB 
PS. 
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equal to the primordial PS (i.e., with each window func- 
tion comprising a single bin). 

We obtain this window matrix for the marginalised 
Fisher matrix of the bins and hence it can only be applied 
to the Fisher matrices of Planck and the combination of 
Planck and SDSS, which are invertiblc. 

8. RESULTS 
8.1. PCA 

First we show the results for the PCA. Principal com- 
ponents, the XiS, obtained for SDSS, Planck and their 
combination are shown as colour-coded matrix plots; X^s 
are shown from left to right with increasing errors (which 

1/2 

is equal to 1/ , as the eigenvalues are constructed for 
the Fisher matrix). Original parameters are shown ver- 
tically starting with the bins on the bottom to #jS on 
the top. For the bins, the vertical width of the box is 
an indication of the binwidth. Some k values for the 
bins are shown on the figures. We group the components 
depending on the Xi, showing those with values greater 
than 0.4; those between 0.2 and 0.4; and those below 0.2. 
SDSS 

The result is shown in Figure [Tj There are a total of 8 
bins that could be obtained to meet the S/N criteria as 
explained above. Together with the 6 0jS, we have a total 
of 14 original parameters and 14 principal components, 
XiS. The last 6 principal components are not measured 
well (they have large/negative eigenvalues — Table [TJ. 
This is because, as explained before, the SDSS Fisher 
matrix is not a positive definite matrix; we have a total 
of 8 data points and this means only 8 parameters (or 8 
different combinations of the parameters, i.e. XiS) can 
be measured. 

The best measured principal component, Xx, has only 
the cosmological parameters (6j) contributing signifi- 
cantly, with h being dominant. The fact that there 
is more than one cosmological parameter contributing 
to this principal component means that SDSS can only 
measure a linear combination of them — a degeneracy 
between these parameters. Xi measures a combination 
of the bins and #;S. Other principal components, X%- 
Xg, measure the bins only, with no contribution from 0,s 
at all, and the correlation amongst the bins is between 
neighbouring ones only. Intuitively, you would expect 
more correlation between the bins and 0<s. Remember 
that the errors for the bins are related to the matter 
transfer function — equation 1101 Therefore, you would 
expect that a change in 9iS would induce a change in 
the matter transfer function and hence a correlation be- 
tween bins and #jS. However, look at Figure [2] where 
it is showing all the derivatives that goes in the Jaco- 
bian. The derivatives with respect to 9iS seem to scale 
relatively close (apart from r where it had to be multi- 
plied by 200). However, the derivative with respect to 
the primordial PS bins has to be rescaled by 10 -8 to fit 
in the same range as the rest of the derivatives. This 
suggests that perhaps the changes in 0iS are not large 
enough in this setting to have a significant effect on the 
matter transfer function and therefore the correlation is 
not that significant to show effects in the PCA. Note that 
the correlation between the bins shows our limits to what 
we can learn about the primordial PS. This correlation 
arises due to our lack of knowledge of the cosmological 
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M.I 10.91 10.81 10.61 10.41 10.21 0.0 

Fig. 1. — The principal components for SDSS with no pri- 
ors on OiS. Xi are shown from left to right with increasing 
errors (= 1/y/Xl). Original parameters are shown vertically 
starting with the bins on the bottom to 8iS on the top. For 
the bins, the vertical width of the box is an indication of the 
binwidth. We group the components depending on the Xi, 
showing those with values greater than 0.4; those between 0.2 
and 0.4; and those below 0.2. (refer to Figure ?? for the colour 
coding). The last 6 principal components can be ignored as 
they are not measured — refer to Table [Tj and text for more 
details. At the bottome we show the colour plot indicating 
different levels of contribution to the principal components. 

parameters. If we knew the parameters perfectly, we 
would have what is shown in Figure which is in fact 
the Fisher matrix itself. Generally it seems that SDSS 
measures cosmological parameters better than the pri- 
mordial PS and within primordial PS bins, it measures 
small scales better than large scales. 

We also investigated what improvements we would see 
given better — realistic — knowledge of the cosmologi- 
cal p arameters. Hence, WMAP5 priors (jDunklev et alJ 
[2008h were added to constrain the 0jS in the Fisher ma- 
trix, by adding the inverse variance of each parameter 
to the Fisher matrix, i.e. ignoring the correlations. The 
result is shown in Figure [H Some of the degeneracies 
between the cosmological parameters have been broken. 
For example, Xi now measures n s almost perfectly Also, 
flj, and r dominate completely in X\\ and Xn respec- 
tively, with no contribution from any other parameter. 
Note that the errors on the principal components have 
reduced and now all X;S, except X\&, can be measured 
well — Table[Tj This is expected as WMAP5 does a good 
job measuring these cosmological parameters. With re- 
spect to bins, it seems that adding priors and improving 
constraints on cosmological parameters has only helped 
to measure linear combinations of the bins better and has 
not been able to break the degeneracy between them. 
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SDSS with priors 



0.00 0.02 0.04 0.06 

k/hMpcT' 

Fig. 2. — The derivative of galaxy PS with respect to the pri- 
mordial PS bins and 6iS, weighted by the parameters values. 
This is exactly what goes in the Jacobian. It is interesting to 
see in this fc-range, much of the variation is on large scales. 
Hence it is no surprise that SDSS measures small scales bet- 
ter. 
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Fig. 3. — The principal components of SDSS for the pri- 
mordial PS bins only, assuming (9;s are known perfectly. No 
correlation exists between the bins. Compare to Figures [T] 
and |4] and see how the lack of knowledge of the cosmological 
parameters induce correlation between (neighbouring) bins. 

Planck 

For Planck there are a total of 23 bins and this, with 
the 6 9iS, means we have 29 principal components, shown 
in Figurc[5j They all seem to be measured well and better 
than SDSS — Table [TJ The reflection of the acoustic 
peaks of CVs on the bin sizes can clearly be seen; the 
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Fig. 4.— The principal components for SDSS with WMAP5 
priors. Now all XiS, apart from X14, are measured well — 
Table[T] Also, some of the degeneracies between 8iS have been 
broken. 

ones corresponding to the peaks are measured with a 
better resolution. To see this, look at equation 



(19) 

which is (S/N) 2 a n. This is equation [5] weighted by the 
signals. First bracket can be ignored as it is almost a 
constant due to the relation F oc C ; ~ 2 — equation[U The 
primordial PS, A^(ki), in the second and third brackets 
can also be ignored as it is a constant. Therefore we are 
left with 

E?^- (20) 



d9 a 89(3 



For the bins the derivative in this equation is the radi- 
ation transfer function as shown in equation 1121 The 
summation over i then gives the oscillatory feature seen 
in k space — see Figure [5] to see the pictorial version of 
this. 

Just like SDSS, Planck seems to measure the cosmolog- 
ical parameters better than the primordial PS and overall 
does a better job than SDSS, giving smaller errors and 
less correlation between them. This is not surprising as 
we already know Planck does a good job measuring the 
cosmological parameters; it measures J1a, h and n s very 
well, with only slight correlation with other cosmological 
parameters. Note that n s is measured almost perfectly 
with no correlation with #,;S (or the bins, as expected). 

The rest of principal components contain the highly- 
correlated bins only, with no particular large contribution 
from any of them. Intuitively one might expect the corre- 
lation to be between neighbouring bins only. The reason 
for the longer-range correlation lies in the form of the 
radiation transfer function; for each i, this transfer func- 
tion spans a fc-range around k= 1/<1a, where d,A is the 
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Fig. 5. — The principal components of Planck. Again bins 
are shown on the bottom and cosmological parameters on the 
top, ordered in the same way shown in previous figures. All 
principal components seem to be measured well and better 
than the SDSS case. Although, no particular scale seem to 
dominate strongly to any of the principal components. 
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Fig. 7. — The derivative of radiation PS with respect to 8iS. 
Note how different parameters dominate on different scales. 
For example, n a dominates on large scales and Qb dominates 
on small scales. If Planck can only measure a linear combina- 
tion of n s and fib, the degeneracy between these parameters 
could induce a degeneracy between large and small scale bins! 
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0.02 0.04 0.06 0.08 0.10 

Fig. 6. — The derivative of r adi ation PS with respect to the 
primordial PS bins, equation 1121 weighted by the primordial 
PS. Note that the bin with I = 400 dominating, gets contribu- 
tions from all Is from 100 to 500. This makes the correlation 
between the bins on all scales possible. 

angular-diameter distance to the last-scattering surface. 
This is due to the projection of a 3D Universe onto a 2D 
sphere around us. Equation [T2l shows what exactly con- 
tributes to the Jacobian for the Fisher matrix analysis. 
For each £, this derivative integrates the radiation trans- 
fer function over the fc-range of the bins. This would be 



reflected as correlation between neighbouring bins. How- 
ever, remember that in the Fisher matrix analysis the is 
get summed over (equation [5]) and this now makes cor- 
relation between all bins possible; Figure [6] shows a pic- 
torial version of equation 1 121 weighted by the primordial 
PS. Note how each £ spans a range of k. The summation 
over all £s means that, for example, the bin with £ = 400 
dominating has contributions from all £s from 100 to 500, 
with each spanning a different range of k. This induces 
correlation between bins of all scales. 

This sort of correlation between small and large scales 
might even be worse when there is a degeneracy between 
the measured cosmological parameters. For example, 
consider an experiment that could only measure a lin- 
ear combination of n s and where n s is dominant on 
large scales and f2t> is dominant on small scales — Fig- 
ure [7] The degeneracy between these parameters could 
induce a degeneracy between large and small scale bins. 

We also investigated if the lack of prior knowledge of 
0jS induces extra correlation between the bins, as in the 
SDSS case. Figure [8] shows the principal components for 
the bins with no OiS — i.e. assuming cosmological param- 
eters are known perfectly. Since Planck's measurements 
of the parameters will be much better than even those 
from WMAP (inclusion of which was able to remove the 
correlations for SDSS), we might expect little change. 
Indeed, not much is changed. The only improvement is 
on the range of errors, which now span a smaller range 
- Table [TJ Note, however, that the smallest error for 
this set is still larger than the smallest error for the set 
including 0jS. This is because 8iS are generally measured 
better than the primordial PS bins and hence they lower 
the errors. Instead, comparing the largest errors of both 
sets shows the improvements. Despite the smaller errors 
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Fig. 8. — The principal components of Planck for the primor- 
dial PS bins only, assuming (9 4 s are known perfectly. Compare 
to Figure [5] It seems like lack of knowledge of cosmological 
parameters does not have much of an effect in measuring the 
primordial PS bins in this case. This could be due to the fact 
that Planck measures 8iS very well. 

for this set, not much is improved in terms of correlation 
between the bins. 

We also want to consider the bins on their own. Hence, 
we consider the correlation between the bins for the 
marginalised Fisher matrix of bins (that is, marginalised 
over the other cosmological parameters, 6*;s). This is ob- 
tained by inverting the parent Fisher matrix to get a 
covariance matrix, which holds the marginalised errors 
for all the parameters. Take the sub-block of this ma- 
trix which holds the errors for the bins and invert this 
to get a marginalised Fisher matrix and diagonalise this 
matrix. The principal components for this Fisher matrix 
are shown in Figure [9l The first thing to note is that 
bins contribute more significantly to some of the princi- 
pal components. In particular there are some mid-scale 
bins which seem to be measured well. For example, look 
at Xig and X22', they seem to have uncorrelated some 
mid-scale bins from the rest of the bins. 

Another interesting result is that very large and very 
small scales never really dominate in the principal com- 
ponents with large errors. They only contribute to them 
at levels of < 0.01. Remember that XiS with large er- 
rors carry the most correlation and therefore the fact 
that mid-scale bins do not contribute to these principal 
components means that they are measured quite well. 

To sum up, it seems like Planck will largely decorrclatc 
the primordial PS from the 6* 2 ;S (and therefore the transfer 
function) but cannot exactly uncorrelate the bins them- 
selves 

Planck fc SPSS 

The results are shown in Figure 1101 Combining sur- 
veys has clearly helped to improve the resolution of the 
primordial PS. Now there are a total of 48 bins in the 
same fc-range. Again the cosmological parameters are 
measured better than the primordial PS and there is 
also less correlation between the cosmological parame- 



Fig. 9. — The principal components of the 'marginalised' 
Fisher matrix of Planck. 
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Fig. 10. — The principal components of Planck & SDSS. 
Clearly, resolution of the primordial PS has improved. Also, 
an almost diagonal trend can be seen now, showing small 
scales are measured better than the large scales. There is 
also less correlation between OiS. 

ters compared to the previous cases. There is also less 
correlation between the bins themselves. Both features 
of SDSS and Planck can clearly be seen here. For exam- 
ple, acoustic oscillations in the C^s still influence the bin 
sizes and resolution of the primordial PS. It also seems 
like small scales are measured better than large scales, 
which is a feature seen in the SDSS case. 

Figure [TT1 shows the results for the marginalised Fisher 
matrix of the bins for SDSS and Planck combined. Com- 
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Fig. 11. — The principal components of the 'mariginalised' 
Fisher matrix of Planck & SDSS. 

pare to Figure [TU1 not much change can be seen. 

8.2. Hermitian Square Root of Fisher Matrix 

Figure shows the window functions for Planck de- 
rived from the Hermitian square root decorrelation. Note 
that only the magnitude of the components of H m s are 
important and not their sign. However, it is worth 
mentioning that for the non-marginalised Fisher matrix 
(both for Planck and its combination with SDSS), these 
window functions have only positive values. Therefore, 
it is the lack of knowledge of the cosmological parame- 
ters (and the induced correlation between the bins) that 
introduces non-physical negative values into the window 
functions. The window functions, H m s, are plotted in the 
order of increasing errors, so that H\ is the best measured 
and H23 the worst measured vector, respectively. Here, 
small scales seem to be measured best and large scales 
measured worst, contributing to H m s with the lowest and 
highest errors respectively. It seems like Planck has not 
been able to decorrelate the bins completely and some 
correlations between neighbouring bins can be seen. In 
addition, bins in the range of k ~ 0.02-0.04/iMpc^ 1 have 
a large contribution to their H m s, compared to the other 
bins. Compare this to Figure fT3l where we diagonaliscd 
the marginalised Fisher matrix through its eigenvectors 
(This is exactly Figure [9] plotted in this format for eas- 
ier comparison). In the PC A case, the correlations seem 
not to be only between neighbouring bins, but between 
bins of all scales, which is not seen in this case! Also, 
the compactness seen here (i.e. more of a window-type 
feature) is not seen in the PCA case; there is no partic- 
ular scale that contributes significantly to the principal 
components. 

Figure [14] shows the windowed PS for Planck. It is 
plotted so that each P m is placed at the k n from which it 
receives the largest contribution. The vertical errors bars 
shown arc A2(fc l )(HF~ 1 H T ), where A? (A*) is the ampli- 



Fig. 12. — The row vectors of H_ for the marginalised Fisher 
matrix of Planck. These vectors are ordered with increasing 
errors, so that Hi is the best and H23 is the worst measured 
vector. This, unlike the principal components, shows that 
correlation is only between neighbouring bins and, that bins 
on large scales are measured better than the ones on small 
scales. 
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Fig. 13. — This is Fig ure [9l plotted in this way for easier 
comparison with Figure 1121 

tude of the primordial PS in the bins and (HF~ 1 H T ) is 
the errors propagated through the H m distribution. The 
horizontal error bars are the half- width at half-maximum 
in each direction of the main peak of each H m . The orig- 
inal primordial PS is plotted for comparison. P m seems 
to be at a lower level than the unwindowed primordial 
PS. Remember that P m is not a physical PS per se. How- 




Fig. 14. — The windowed PS obtained from Planck. 

ever, the observed differences from the original PS arise 
due to the induced correlations between the bins. In a 
perfect setting, where there are no correlations between 
bins, you do expect P m = A^(fc) to be true. Note that 
the main feature of this plot is that vertical errors, unlike 
those for the original primordial PS, are not correlated. 
The correlation between the errors has been transferred 
to overlaps between the window functions — as shown 
in Figure [j"2l 

Figures [TBI and [T71 show the same set of results for com- 
bination of Planck and SDSS. Again large scales are con- 
tributing to H m s with the largest errors. There seems to 
be less correlation between neighbouring bins compared 
to the Planck case. Also, note that bins in this case are 
narrower and therefore correlation between neighbouring 
bins still means correlation between a narrower range 
of k. Compare Figure [15] to Figure [16] (same as Fig- 
ure [TTj) . Again, there is less compactness in the PCA 
case, however more than what is seen for Planck on its 
own. Figure [15] indcates that bins in the vicinity of 
k ~ 0.02 — 0.025/iMpc _1 seem to contribute very strongly 
to their H m s compared to other bins, in particular the 
last window function, H^. This effect gets carried on to 
P rn , with Pn having a very large amplitude — Figure 

EH 

9. CONCLUSIONS 

The primordial PS holds precious information about 
the physics of the early Universe and constraining it has 
been one of the key goals of the modern cosmology. How- 
ever, the induced degeneracy between the cosmological 
parameters of the matter/radiation transfer functions 
and the primordial PS limit our ability to recover the 
primordial PS completely, even from a perfect survey, 
espec ially in the case of CMB surveys (]Hu fc Okamotol 
2004). Different surveys probe different scales with dif- 
ferent accuracies and might not be able to constrain 
the primordial PS to a desired resolution on their own. 



Fig. 15. — The row vectors of H_ for the marginalised Fisher 
matrix of Planck & SDSS. Like before, they are ordered with 
increasing errors. The correlation between neighbouring bins 
still exists but to a lesser extent. Also, note that the bins are 
narrower here so that correlation between neighbouring bins 
still means a correlation within a narrower fc-range. 

Planck & SDSS 
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Fig. 16. — This is Fig ure [TT1 plotted in this way for easier 
comparison with Figure 1151 

However, put together, they make significant improve- 
ments. In this paper we have investigated these lim- 
its/improvements for Planck and SDSS. For this purpose, 
we have assumed a non-parametric function of the pri- 
mordial PS and have constructed a parameter space con- 
taining a set of carefully chosen bins of the primordial 
PS along with a set of cosmological parameters. We con- 
structed a Fisher matrix for this parameter space for the 
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Fig. 17. — The windowed PS obtained from combination of 
Planck & SDSS. 

two different surveys separately and combined. By diago- 
nalising these Fisher matrices, via two different methods 
of eigenvector decomposition (PCA) and the Hcrmitian 
square root, we have investigated the induced correlation 
between the primordial PS bins and the cosmological pa- 
rameters. 

In the PCA case, we came to conclude that SDSS and 
Planck together measure the cosmological parameters to 
a better extent, and even break the degeneracy between 
them. They can increase the resolution of the primordial 
PS by about twice as much and can also condense the 
correlation between bins to be only amongst neighbour- 
ing ones. On the whole it seems like they can constrain 
small scales better than large scales. 

By the use of Hcrmitian square root of the Fisher ma- 
trix we managed to divert the correlation amongst the 
marginalised errors of the bins to the correlation between 
the bins themselves. In this case, combination of SDSS 
and Planck helped to decrease the level of correlation be- 
tween neighbouring bins, but also, because it has helped 
to increase the resolution of the bins, correlation between 
neighbouring bins means correlation between a smaller 
range of k. 

Clearly adding the two surveys have helped to con- 
strain the primordial PS to a better degree. Obviously, 
further surveys such as Ly-a (e.g. SDSS LyaF PS), 
weak lensing (e.g. Euclid), peculiar velocity (e.g. Cluster 
Imaging Experiment (CIX)), etc. can help even more to 
measure the primordial PS. 
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TABLE 1 

Errors for different sets for SDSS, Planck and 

COMBINATION OF PLANCK AND SDSS. 



SDSS 


x 1 


x 2 


x 3 


x 5 


X 7 


x 9 


X u 


X 13 


A14 


No priors 


0.0038 


0.0160 


0.0287 


0.0328 


0.0339 


7E5 


NaN 


NaN 


NaN 


No Oi 


0.0282 


0.0299 


0.0317 


0.0340 


0.0359 










WMAP5 priors 


0.0038 


0.0123 


0.0236 


0.0308 


0.0340 


0.0357 


0.0709 


0.02571 


23.38 



Planck 


x 1 


x 2 


x 3 


x 5 


Xio 


X15 


X20 


x 25 


A28/29 


PCA-No priors 
PCA-No fts 
PCA-Margin. 
Hermitian Sqrt 


0.0004 
0.0110 
0.0204 
0.0236 


0.0022 
0.0149 
0.0210 
0.0261 


0.0035 
0.0181 
0.0255 
0.0268 


0.0152 
0.0224 
0.0313 
0.0283 


0.0295 
0.0333 
0.0410 
0.0386 


0.0402 
0.0422 
0.0553 
0.0487 


0.0515 
0.0572 
0.0890 
0.2525 


0.0700 


0.4953 



Planck & SDSS 


Xi 


X 2 


x 3 


Xio 


X20 


X30 


X40 


X50 


X54 


PCA-No priors 
PCA-No 9 iS 
PCA-Margin. 
Hermitian Sqrt 


0.0004 
0.0561 
0.0254 
0.0323 


0.0020 
0.0568 
0.0289 
0.0325 


0.0035 
0.0575 
0.0304 
0.0327 


0.0348 
0.0624 
0.0425 
0.0374 


0.0463 
0.0729 
0.0523 
0.0578 


0.0548 
0.0891 
0.0598 
0.0762 


0.0642 
0.1003 
0.0823 
0.2102 


0.1132 


0.5261 
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